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ABSTRACT 



A light pseudoscalar coupled to two photons would be copiously emitted by the 
core of a supernova. Part of this flux would be converted to 7— rays by the galactic 
magnetic field. Measurements on the SN1987A 7— ray flux by the Gamma-Ray 
Spectrometer on the Solar Maximum Mission satellite already imply a bound on the 
coupling g < 3 x 10 -12 GeV -1 . The improved generation of satellite-borne detectors, 
like EGRET or the project GLAST, could be able to detect a pseudoscalar-to- 
photon signal from a nearby supernova, for allowed values of g. 
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Introduction 



Pseudoscalar particles are often fundamental ingredients of Particle Physics models. Examples 
are axions |I[ or Majorons 0], coming from models with spontaneous breaking of a Peccei- 
Quinn symmetry || or of a global lepton symmetry, respectively. Other examples are light 
bosons from extra-dimensional gauge theories and arions || . Pseudoscalar particles usually 
couple to two photons via the interaction lagrangian 

C=\g4> e^pF^F^. (1) 

Limits on the coupling g come from laboratory experiments and from cosmological and 
astrophysical observations. These have been discussed and collected in ||. Further astro- 
physical constraints have been recently examined by Mori |7j]. The limits on g depend on the 
mass m of the pseudoscalar. For very light masses, m < 10 -9 eV, the best constraints come 
from astrophysics (such constraints are not interesting for axion models, where g and m are 
related; but are of interest to other models where m < 10~ 9 eV does not imply an exceedingly 
small coupling g). Mohanty and Nayak || considered the creation of a background due to 
strong magnetic fields in pulsars. Pulsar signals propagating through this background show 
a time lag between different modes of polarization. They found the limit 

g < 5.3 x 10- u GeV -1 , (2) 

valid for m < 10~ 10 eV. Another way to constrain g for m < eV has been discussed by 
Carlson || . He studied x-ray conversion of <fi produced in giant cores and found the following 
limit 

9 <2.5x KT 11 GeV~\ (3) 
using HEAOl satellite data on a-Ori x-ray emission. A stringent bound on g has been recently 



obtained by Krasnikov [10]: conversion would lead to a large scale anisotropy of the 



x-ray background. The observed isotropy leads to the estimation [10 



g<lx GeV-\ (4) 

This author also considers similar effects to those studied by Carlson ||. In addition he dis- 
cusses polarimetric effects in the emission by magnetic white dwarfs, and claims measurements 
may be sensitive to g ~ 10~ n GeV~ n . 

In the present paper, we show that it is possible to improve these bounds on g for m < 10~ 9 
eV. We will use the fact that, should these pseudoscalars exist, they are copiously produced 
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when a supernova occurs. For small g, these particles will stream away from the supernova 
core, without further interactions. In their path, the galactic magnetic field will convert a 
fraction of the flux back to gamma-ray photons. At the time the neutrino burst from SN1987A 
was observed, the Gamma-Ray Spectrometer (GRS) on the Solar Maximum Mission (SMM) 
satellite was operative and did not observe such a gamma-ray signal. We use this null result 
to set the stringent limit 

g<3x KT^GeV -1 . (5) 

The next generation of detectors, like EGRET on the Compton GRO satellite or the project 
GLAST, could be able to detect a 7— ray signal due to <ft conversion from a nearby supernova 
collapse, for values of g allowed by our analysis. 



Pseudoscalar production 

Immediately after collapse, <fi is produced in the hot and superdense core of the supernova. In 
what follows, we take a core temperature T = 60 MeV, proton number density n p = 1.4 x 10 -38 
cm -3 and radius R c = 10 km. As has been discussed in ||11|| , a source of uncertainty arises 



because we do not know which is the equation of state at supernuclear densities. To estimate 
the uncertainty, we will follow that reference and change the central density by a factor of 
2. All the other parameters will change correspondingly. For example, the core temperature 
changes by a factor 1.6, if one assumes an isentropic collapse. Notice that this allows the core 
temperature to span the conventional range of T ~ 30 MeV — 100 MeV. Since couples to the 
electromagnetic field, the relevant interactions will involve electrons and protons. The Fermi 
momentum of electrons and protons is pp = 320 MeV. For the relativistic electrons, Pf/T 3> 1 , 
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thus they are extremely degenerate. As to the non-relativistic protons, ^-/^p — 0.6, and 
hence they are only moderately degenerate. 

Let us start considering the Primakoff process on protons, py — ► pcj), that leads to <fi 
creation. The number of pseudoscalars produced per unit volume and per unit time, at 
temperature T and with energy between E min and E max is 



- ma :r 



N(T) = / a(u)v n p dn 7 (T, u). (6) 



En 



We integrate between E min and E max since gamma-ray detectors are only sensitive to a fixed 
energy band. In expression (^|) n p and n 7 are the number densities of protons and photons 
respectively, v their relative velocity and the Primakoff cross section as a function of the 
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photon energy u is (see [O 



a{u) = \ag 2 



W 1 + K -i 



4u; 2 / \ k 



2 



(7) 



The screening wave number k appears because of the collective behavior of the plasma, which 
cuts off the range of the Coulomb potential for scales larger than ~ re -1 . It is given by 

K — Kp -\- Krpp (8) 

4 = ^ (9) 

K TF = PfE F , (10) 

7T 

with ^op and Ep the Fermi momentum and energy of the electrons. The proton contribution 
in our case is k,d = 50 MeV, where we are considering the protons not degenerate (see below). 
The contribution of the degenerate electrons is given by ktf = 33 MeV, smaller than the 
proton contribution. Protons can move in the plasma (or better said, in momentum space) 
more freely than the degenerate electrons and so they can screen charges more easily than 
the electrons. 



In order to obtain a(u) in eq. (|7|) one needs to know the thermal average < \F(q)\ 2 >, 
being F(q) the usual form factor associated to the charge distribution of the plasma. To make 
the thermal average we need pfj(r), the probability per unit volume of finding the charge Zj 
at a distance r of the charge Z{ [O 



(u: 



for an arbitrarily large volume V. In addition to this charge correlation one also has to consider 
the statistical correlation between fermions. The probability pfj{r) of finding a fermion j at 
a distance r of the fermion i is [[y| 

P?j(r) = ^ (l - exp(-mTr 2 )) . (12) 

One may estimate the relative importance of the two different types of correlation by evalu- 
ating the quotient 

V~ l — pfAr) 

e = jk at r = k . (13) 

V-i-pg(r) 

Using the physical parameters in the supernova core we obtain e < 10 -5 for the gas of 
protons. Therefore we can neglect the near degeneracy of protons and consider only the 
charge correlation. In this case, one gets expression (01). 



One might also consider other mechanisms of pseudoscalar production such as the Pri- 
makoff process on electrons, — ► e<p, or processes involving only nucleons in the initial 
state, pn — > pwycj), where a virtual photon attached either to a proton or to a virtual 
charged pion splits into the final 7 and <fi. These processes would add to the one previously 
considered and would make our limit on g more stringent. However, the process e7 — > e<p 
is suppressed compared to py — > p<p due to the extreme degeneracy of the electrons; pro- 
duction of a pseudoscalar entails a change of momentum Ap e ~ T which is not allowed for 
the bulk of electrons in the Fermi sea (one has to bear in mind that forward peaks, in which 
Ap e ^ 0, are suppressed by the finite range of the Coulomb potential in a plasma). Processes 
like pn — > pwy(f) are also unimportant since they can be visualized as 70 production by the 
nearly static electric field created by protons. This sort of production is clearly diminished 
by energy conservation. 

Finally, we can write the expression for the pseudoscalar flux on the Earth, at a distance 
D = 55 kpc from SN1987A, considering pj — > p<p as the only production mechanism in the 
supernova core: 

g \ 2 / 55kpc \ 2 / gg_W Hp )( T 
10- n GeV~V I D ) UOkm/ U.4 x 10 38 cm- 3 / UoMeV, 



3 x 10 4 cm-V 1 



(14) 



where 



and being 



y(C ) ^mirii E max ) 



2tt 



dx- 



e x -1 



y + C) ln(l + x 2 /C) - x 



k/2T, 



•^min ^ rain l^i 



(15) 

(16) 
(17) 
(18) 



One can check that the energy drain by these pseudoscalars during the collapse of the 
supernova core is at least a factor 10~ 4 smaller than the energy released in neutrinos, therefore 
we can ignore the back reaction caused by <fi emission on the supernova evolution. 



Pseudoscalar conversion 



The mixing between the photon and low mass particles in magnetic fields leads to very inter- 
esting phenomena. Since the pioneering work in jTJJ], a variety of implications for laboratory 
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experiments and astrophysical observations have been investigated (see references in [El). 



For our purposes, we use the formalism developed in [15]. The <p ~ > 7 transition prob- 
ability, for a beam traversing a transverse magnetic field B? after a distance L, is given 
by 



P{^l) = \fBll? {^-)\ (19) 



x 



with 



and 



x = -^Al c + g*B 2 T (20) 



^osc 



Up — m 2 



2uo 



(21) 



where a; is the energy and 
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2 - We - (6.4 x ID- 2 eV) 2 (—^) (22) 



p m e \ ' VO.03 cm 

is the plasmon mass. We have normalized it to the mean electron density in the interstellar 
medium, n e ~ 0.03 cm 3 . 



A coherent effect, which implies sina;/a; — ► 1, is obtained provided 

m < 10" 9 eV, (23) 

so that our constraint on g will be valid only for such small masses. In order to evaluate the 
probability in (0), we need to specify the magnetic field structure. We will adopt the model 
used in 0, consisting of a toroidal 2/xG magnetic field. The coherence length is about the 
order of several kiloparsecs jnS |. To be conservative we take a coherence length L = 1 kpc. 



In the direction of SN1987A one has 

1 /2 

B T = 2/iG (l - sin 2 I cos 2 b) ~ 1 //G, (24) 

where the corresponding galactic coordinates have been used. 
Summing up, we have 



10- u GeV"V \lfiGl \lkpc 



Results and Conclusions 



At the time the neutrino burst from SN1987A reached the Earth, the satellite-b orne GRS 
was on duty measuring the incident 7— ray flux. This measurement provides an observational 
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limit on the 7— ray flux coming from the supernova, and consequently on the photons from 
supernova emission. Indeed, data from SMM [|l7] in the energy band E m i n = 25 MeV 
< E < E max = 100 MeV imply that 

$ P((f> -»• 7) At < 0.6 cm" 2 . (26) 



The characteristic time over which the proto-neutron star emits most of its gravitational 
energy and therefore the bulk of the hypothetical (f) particles is the diffusion time of neutrinos. 
Detailed stellar evolution calculations |1| render about 10 seconds for this diffusion time scale. 
To be conservative we shall take At = 5 s which is roughly the characteristic signal decay 
time for the neutrino burst of SN1987A. We thus find the limit 

g<3x 10- 12 GeV" 1 . (27) 

Changing the physical parameters of the supernova as previously discussed, we estimate an 
uncertainty factor of 2 in the limit 



Finally, we would like to make some comments on the energies of the photons from pseu- 
doscalar conversion. The expected spectrum has a mean energy value slightly bigger than 
2.7T ~ 160 MeV (the energy spectrum is not exactly that of a black body since the Pri- 
makoff cross section increases with the energy). Thus, it extends above the GRS cut-off at 
E m ax = 100 MeV and, consequently, part of it would not have been detected. Better prospects 
could be expected with EGRET on the Compton GRO satellite, launched in 1991, since it is 
able to detect 7— rays from E min = 20 MeV up to E max = 30 GeV and therefore is sensitive 
to the whole spectrum. In addition, the EGRET detector is more sensitive than the GRS 



detector. A pseudoscalar with g allowed by equation (p7|) would possibly give a 7— ray sig- 
nal in EGRET from a nearby supernova. Even more promising is the project GLAST with 
detectors to measure 7— rays in the energy range from E min = 50 MeV to E max = 100 GeV, 
and a factor 100 better in sensitivity than EGRET. 
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